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I.  INTRODUCTION 


One  of  the  most  important  aerodynamic  performance  characteri sti cs  for 
shell  is  the  total  drag.  The  total  drag  for  projectiles  can  be  divided  into 
three  components:  (1)  pressure  drag  (excluding  the  base),  ( i i )  viscous  (skin 
friction)  drag,  and  (i 1 i )  base  drag.  The  base  drag  is  a  major  contributor  to 
the  total  drag,  particularly  at  transonic  speeds.  Thus,  the  determination  of 
base  pressure  is  essential  in  predicting  the  total  drag  for  projectiles.  The 
breakdown  of  the  total  drag  Into  various  components  is  important  in  the  pre¬ 
liminary  design  stage  of  shell.  This  information  can  aid  the  designer  to  find 
potential  areas  for  drag  reduction  and  achieve  a  desired  Increase  in  range 
and/or  terminal  velocity  of  projectiles. 

Design  codes  are  available  that  can  predict  the  individual  drag  compo¬ 
nents  and  thus,  the  total  aerodynamic  drag.  Typically,  the  total  drag 
predicted  by  these  design  codes  agrees  well  with  the  total  drag  measured  from 
flight  tests.  Rut,  how  well  do  these  codes  predict  the  individual  drag  compo¬ 
nents  for  a  projectile?  It  is  difficult  to  measure  the  individual  drag  compo¬ 
nents  experimentally  (especially  the  base  drag  and  skin  friction  drag)  and 
verify  the  accuracy  of  these  predictions.  Recently  developed  Navi er-Stokes 
computational  procedures  are  capable  of  predicting  all  the  individual  drag 
components  for  projectiles  and  can  be  used  to  determine  how  well  the  design 
codes  predict  the  individual  drag  components.  This  report  describes  such  a 
Navi er-Stokes  computational  study  to  predict  the  individual  drag  components  of 
projectiles  and  determine  the  accuracy  of  the  predictions  from  available 
design  codes. 


II.  COMPUTATIONAL  TECHNIQUE 

The  Azimuthal  Invariant  (or  Generalized  Axi symmetri c)  thin-layer  Navier- 
Stokes  equations  for  general  spatial  coordinates  S,  n,  4  can  be  written  as1 

a  q  +  3  E  +  3.G  +  H  =  Re"  3  S  (1) 


where 


and 


5  =  4(x,y,z,t)  is  the  longitudinal  coordinate 
n  =  n(y,z,t)  is  the  ci rcumferenti al  coordinate 
C  =  c(x,y,z,t)  is  the  near  normal  coordinate 
t  *  t  is  the  time 


2»  C.J.  Nietubicz ,  T.H.  Pulliam ,  and  J.L.  Steger,  "Numerical  Solution  of  the 
Azimuthal-Invariant  Uavier-Stokes  Equations US  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  AR3RL-TP-0222?  ,  March 
1980.  (AD  A085718)  (Also  see  AIAA  Journal,  Vol.  13,  Ro.  12,  December 
1980  ,  pp.  1411-14 12 1 
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U  *  +  ^xu  +  5yv  +  Czw 

V  =  nt  +  nxu  +  nyv  +  n2w  (2) 

W  =  C.  +  C  u  +C  v  +  c  w 
t  x  y  z 

represent  the  contravari ant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  nondimensi onal i zed  with 
respect  to  aa  (free  stream  speed  of  sound).  The  density  (p)  is  referenced 
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to  and  total  energy  (e)  to  p^a2.  The  local  pressure  is  determined  using 
the  equation  of  state, 

P  =  (y  -  1 ) [e  -  0.5p(u2  +  v2  +  w2)]  (3) 


where  y  is  the  ratio  of  specific  heats. 

In  Equation  (1),  axisymmetric  flow  assumptions  have  been  made  which 

result  in  the  source  term,  H.  The  details  of  how  this  is  obtained  can  be 
found  in  Reference  1  and  are  not  discussed  he^c.  Equation  (1)  contains  only 
two  spatial  derivatives.  However,  it  ret'ins  all  three  momentum  equations  and 
allows  a  degree  of  generality  over  the  standard  axisymmetric  equations.  In 
particular,  the  ci rcumferenti al  velocity  is  rv  assumed  to  be  zero,  thus 
allowing  computations  for  spinning  projectiles  to  he  accomplished. 

The  numerical  algorithm  used  is  the  Beam-Warming  fully  implicit,  approxi¬ 
mately  factored  finite  difference  scheme.  The  algorithm  can  be  first  or 
second  order  accurate  in  time  and  second  or  fourth  order  accurate  in  space. 
Since  the  interest  is  only  in  the  steady-state  solution,  Equation  (1)  is 
solved  in  a  time  asymptotic  fashion  and  first  order  accurate  time  differencing 
is  used.  The  spatial  accuracy  is  fourth  order.  Details  of  the  algorithm  are 
Included  in  References  2-4. 

For  the  computation  of  turbulent  flows,  a  turbulence  model  must  be  sup¬ 
plied.  In  the  present  calculations,  the  two  layer  algebraic  eddy  viscosity 
model  developed  by  Baldwin  and  Lomax5  is  used.  In  their  two-layer  model,  the 
inner  region  follows  the  Prandtl-Van  Driest  formulation.  Their  outer 
formulation  can  be  used  in  wakes  as  well  as  in  attached  and  separated  boundary 
layers.  In  both  the  inner  and  outer  formulations,  the  distribution  of 
vorticlty  is  used  to  determine  length  scales,  thereby  avoiding  the  necessity 
of  finding  the  outer  edge  of  the  boundary  layer  (or  wake).  The  magnitude  of 
the  local  vorticity  for  the  axisymmetric  formulation  is  given  by: 


2.  J.L.  Stager,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoils AIAA  Journal,  Vol.  IS, 
Ho.  7 ,  July  1978,  pp.  679-686. 

3.  T.H.  Pulliam  and  J.L.  Steger,  "On  Implicit  Finite-Difference  Simulations 

of  Three-Dimensional  Flow,"  AIAA  Jownxal,  Vol.  18,  Ho.  2,  February  1980 , 
pp.  169-167.  — 

4.  R.  Beam  and  R.F.  Harming,  "An  Implicit  Factored  Scheme  for  the  Compres¬ 

sible  Havier-Stokes  Equations,"  AIAA  Journal,  Vol.  16,  Ho.  4,  April  197  9, 
pp.  393-402.  ~~  '  — — 

5.  B.S.  Baldwin  and  H.  Lomax,  "Thin-Layer  Approximation  and  Algebraic  '!odel 
for  Sepaxxxted  Turbulent  Flows,"  AIAA  Paper  Ho.  78-257 1  1978. 
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In  determining  the  outer  layer  length  scale  a  function5 


F(y)  -  y M  [1  -  exp(-y+/A+)]  (5) 

was  used  where  y+  and  A+  are  boundary  layer  parameters. 

The  thin-layer  Navi er-Stokes  computational  technique  described  above  has 
been  used  in  conjunction  with  a  unique  flow  field  segmentation  proce¬ 

dure6,7  which  allows  the  entire  flow  field  over  a  projectile  including  the 
base  region  flow  to  be  computed  An  important  advantage  of  this  segmentation 
procedure  lies  in  the  preservation  of  the  sharp  corner  at  the  base;  in  other 
words,  no  approximation  or  rounding  of  the  actual  sharp  corner  at  the  base  is 

made.  The  details  of  this  procedure  can  be  found  in  References  6  and  7. 

Since  the  entire  flow  field  over  the  projectile  is  calculated,  all  the 

individual  drag  components  can  be  computed  and  thus,  the  total  aerodynamic 
drag  can  be  determined. 


III.  DESIGN  CODES 


A.  NSWCAP 


The  NSWCAP  code8,9  is  a  semi -empi rical /analyti cal  technique  which 
provides  fast  predictions  of  static  and  dynamic  coefficients  of  shell  at 
transonic  and  low  to  moderate  supersonic  velocities.  This  code  is  relatively 
simple  to  use  and  is  formulated  to  be  a  ‘design  tool1. 


6.  J.  Sahu,  C.J.  Nietubicz ,  and  J.L.  Steger,  "Numerical  Computation  of  Base 
Flou  for  a  Projectile  at  Transonic  Speeds,"  US  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground ,  Maryland ,  AR8RL-TR-02495  ,  June  1983. 
(AD  A130293)  (Also  see  AIAA  Paper  Ho.  82-1338 ,  August  1982) 

7.  J.  Sahu ,  C.J.  Nietubicz ,  and  J.L.  Steger,  "Navier-Stokes  Computations  of 
Projectile  Base  Flou  with  and  without  Base  Injection,"  US  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02532 , 
November  1983.  (AD  A135738 )  (Also  see  AIAA  Journal,  Vol.  23,  No.  9, 
September  1985 ,  pp.  1348-1355) 

8.  F.G.  Moore  and  R.C.  Swanson ,  "Aerodynamics  of  Tactical  Weapons  to  Mach 
'Dumber  3  and  Angle- of- Attack  15° ,  Part  I  -  Theory  and  Application," 
NSWC/DL  TR-3S84  ,  February  1977. 

9.  F.G.  Moore  and  R.C.  Swanson,  "Aerodynamics  of  Tactical  Weapons  to  Mach 
Number  J  and  Angle- of- Attack  15°,  Part  II  -  Computer  Procram  and  Usage," 
NSWC/DL  TR  3600 ,  March  1977. 
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The  total  drag  is  obtained  by  a  linear  superposition  of  pressure  drag, 
skin  friction  drag  and  base  drag.  This  code  also  has  the  ability  to  estimate 
the  drag  due  to  a  rotating  band  if  it  is  present. 

The  pressure  drag  consists  of  contributions  from  the  nose  and  the  boat- 
tail.  The  nose  drag  is  based  on  empirical  as  well  as  experimental  data  at 
transonic  speeds.  The  boattail  pressure  drag  is  based  upon  a  small  distur¬ 
bance  potential  solution.  The  estimation  of  nose  drag  and  boattail  drag  at 
supersonic  speeds  is  based  on  the  Van  Dyke  second-order  theory10. 

The  skin  friction  component  of  the  total  drag  is  computed  using  the  model 
of  Van  Driest11  The  base  drag  prediction  is  empirical.  It  is  assumed  that 
the  boattail  is  located  after  a  relatively  long  afterbody  so  that  the 
approaching  external  flow  is  at  free  stream  conditions.  The  base  drag  is 
determined  by  the  expression 


CAB  ■  -Cp8A("-><VVf>J 

where  CD  (M  )  is  the  base  pressure  coefficient  for  a  long  afterbody  with  no 
RA 

boattail.  Rg  is  the  base  radius  and  Rref  is  the  reference  body  radius 

(usually  the  radius  of  the  cylindrical  section).  CD  is  based  on  the  data 

BA 

for  a  long  cylindrical  afterbody  and  a  fully  turbulent  boundary  layer  ahead  of 
the  base.  The  effect  of  base  bleed  or  rocket  exhaust  can  not  be  accounted  for 
in  this  approach. 

B.  MCDRAG 

This  is  another  program12  based  on  semi -empi ri cal  technique  which 
provides  a  ouick  response  and  is  very  easy  to  use.  This  code  is  used  for 
estimating  the  drag  of  a  projectile  in  the  Mach  number  range  of  0.5  to  5.0. 
The  total  drag  takes  the  form 


10.  M .J .  Van  dykes,  "The  Similarity  Rules  for  Second-Order  Subsonic  and 
Supersonic  Flou,"  RAC A  Tech  Rote  38?  5,  October  1958. 

11.  F.R.  Van  Driest,  "Turbulent  Boundary  Layers  in  Compressible  Fluids," 
Journal  of  the  Aeronautical  Sciences ,  V ol.  18,  Ro.  3,  1951 ,  pp.  145-180. 

12.  R.L .  McCoy,  "McDrag  -  A  Computer  Program  for  estimating  the  Drag 
Coefficient  of  Projectiles US  Army  Ballistic  Research  Laboratory , 
Aberdeen  Proving  Ground,  Maryland,  AR3RL-TR-02293 ,  February  1981.  I  AD 
AO  9 8110) 
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where 


Cn  »  total  drag  coefficient  at  zero  angle  of  attack 


pressure  drag  or  wave  drag  coefficient 


viscous  or  skin  friction  drag  coefficient 


base  drag  coefficient. 


The  pressure  drag  consists  of  drag  due  to  the  projectile  nose,  boattail  and 

the  rotating  band  (If  any).  Prediction  of  the  projectile  nose  drag  Is  based 

on  analytical  theories  at  supersonic  speeds  while  the  transonic  nose  drag  Is 
based  on  correlations  with  experimental  data.  This  Is  similar  to  the  pro¬ 
cedure  used  In  the  NSWCAP  code.  The  effect  of  leading  edge  bluntness  Is 
accounted  for  In  estimating  the  nose  drag.  As  for  the  boattail  drag,  second 
order  theory  Is  used  to  correlate  with  experimental  data  to  estimate  this 
component  of  drag  at  supersonic  speeds  and  a  similar  procedure  Is  used  at 

transonic  speeds  as  well.  The  drag  due  to  the  rotating  band  Is  usually  very 

small  (of  the  order  a  few  percent  of  the  total  drag)  and  Is  based  on 
correlations  with  few  experimental  tests. 

The  skin  friction  drag  coefficient,  Cn  ,  Is  given  by. 


where  Cp  *  skin  friction  coefficient  for  a  smooth  flat  plate 
Sy  »  projectile  wetted  surface  area  (excluding  the  base) 


Cp  Is  estimated  analytically  depending  on  whether  the  boundary  layer  Is 

laminar  or  turbulent  on  the  nose.  The  flow  over  the  afterbody  Is  assumed  to 
be  turbulent  always. 

The  base  drag  coefficient,  Cn  Is  estimated  empirically.  The  approach 

B 

taken  here  differs  from  the  procedure  used  in  the  NSWCAP  code.  The  pressure 
drag  and  skin  friction  drag  are  estimated  as  described  previously.  These  con¬ 
tributions  are  then  subtracted  out  from  the  measured  total  drag  coefficients 
which  are  available  from  free  flight  data  of  various  projectiles.  An  average 
base  pressure  Is  then  Inferred  from  the  derived  base  drag  coefficient.  The 
estimate  of  base  drag  coefficient  Is  obtained  from  the  relation, 
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B  yM„ 


(1  -tt) 
00 


where  Pg/P„  a  base  pressure  obtained  from  a  least  square  fit  of  the  data  which 
includes  Mach  number  and  boattal 1  effects 


and 


dB 


■  diameter  at  the  base. 


IV.  RESULTS 

Computations  of  various  drag  components,  and  thus  the  total  drag  have 

been  made  for  three  projectiles:  SOC;  SOCBT;  and  the  M549.  For  the  Navier- 

Stokes  computations,  solutions  were  marched  In  time  until  the  steady  state 
results  were  achieved.  Results  are  now  presented  for  the  three  cases 
considered. 

A.  SOC  Configuration 

This  projectile  has  a  3  caliber  (1  caliber  =  1  reference  body  diameter) 
secant-ogive  nose  and  a  3  caliber  cylindrical  afterbody  as  shown  In  Figure 
1.  One  of  the  first  steps  before  performing  Navler-Stokes  calculations  Is  the 
determination  of  a  computational  grid.  An  expanded  view  of  the  grid  near  the 
projectile  Is  shown  in  Figure  2.  The  grid  consists  of  114  points  In  the 
longitudinal  direction  (Including  30  points  In  the  base  region)  and  50  points 

In  the  normal  direction.  The  grid  points  In  the  normal  direction  are 

stretched  away  from  the  surface  exponentially.  The  clustering  near  the  body 
surface  Is  required  to  resolve  the  boundary  layer.  Additional  grid  clustering 
Is  used  In  the  longitudinal  direction  near  the  nose-cylinder  junction  and  the 
base  where  appreciable  changes  In  the  flow  field  variables  are  expected. 
Figure  3  shows  an  expanded  view  of  the  grid  in  the  base  region.  Fifty  grid 
points  are  used  In  the  base  region  from  the  base  corner  down  to  the  center 
line  of  symmetry.  This  Is  done  to  provide  adequate  grid  resolution  along  the 
base  of  the  projectile. 

The  individual  drag  components  (pressure  drag,  viscous  drag  and  base 
drag)  are  obtained  from  the  computed  solutions  and  their  variations  with  Mach 
number  are  presented  next.  These  results  are  compared  with  the  predictions 
from  the  design  codes.  As  shown  In  Figure  4,  the  pressure  drag  predicted  by 
the  Navler-Stokes  code  agrees  well  with  the  MCDRAG  and  NSWCAP  predictions 
except  at  M  -  1.2  where  a  small  discrepancy  is  found.  Comparison  of  viscous 
drag  Is  shown  In  Figure  5.  Here  the  viscous  drag  predicted  by  the  Navler- 
Stokes  code  Is  In  very  good  agreement  with  the  MCDRAG  prediction.  The  NSWCAP 
code  underpredicts  this  drag  by  about  15%  in  the  transonic  speed  regime.  This 
contribution  to  the  total  drag  Is,  however,  small. 

Figure  6  shows  the  variation  of  base  drag  as  a  function  of  Mach  number. 
As  expected,  the  base  drag  decreases  as  Mach  number  Is  Increased  from  1.2  to 
3.  For  this  range  of  speeds  the  Navler-Stokes  code  result  agrees  very  well 
with  the  NSWCAP  prediction.  MCDRAG  overpredicts  these  results  by  about  12%. 
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Comparison  of  the  base  drag  Is  shown  more  clearly  In  Figure  7.  The  Navler- 
Stokes  results  are  Indicated  by  circles,  experimental  results13  by  triangles, 
data  base14  results  by  diamonds  and  MCDRAG  results  by  squares.  The  results 
from  the  data  base  are  based  on  correlation  of  base  pressures  obtained  from  a 
number  of  experiments  and  other  analytical  techniques.  The  increase  In  base 
drag  with  Increase  In  Mach  number  In  the  transonic  regime  Is  predicted  by  all 
techniques.  However,  large  scatter  In  the  predictions  of  base  drag  exists. 
The  Navler-Stokes  result  agrees  well  with  the  result  from  the  data  base  except 
at  M  ■  1.1  where  It  overpredicts  the  base  drag  by  125.  The  MCDRAG  prediction 
of  base  drag  Is  consistently  higher  (up  to  105)  than  the  data  base  result. 

The  total  aerodynamic  drag  Is  obtained  from  the  Individual  components 
and  is  plotted  as  a  function  of  Mach  number  In  Figure  8.  The  Navler-Stokes 
result  compares  very  well  with  NSWCAP  prediction  at  all  speeds  (transonic  and 
supersonic).  The  MCDRAG  prediction  Is  in  reasonable  agreement  with  other 
results  at  higher  Mach  n  ibers  (M  >  1);  however,  discrepancy  up  to  20%  is 
found  at  low  Mach  numbers  (.9  <  M  <  1.0).  The  total  drag,  as  expected, 
Increases  as  Mach  number  Increases  from  0.9  to  1.2  In  the  transonic  speed 
regime. 


8.  SOCBT  Configuration 

The  model  geometry  for  this  secant-ogl ve-cyl 1 nder-boattal 1  (SOCBT)  con¬ 
figuration  Is  shown  In  Figure  9.  It  has  a  3  caliber  secant-ogive  nose,  a  2 

caliber  cylinder  and  a  1  caliber,  7°  boattall.  A  computational  grid  was  ob¬ 

tained  for  this  configuration  and  an  expanded  view  of  the  grid  near  the  pro¬ 
jectile  Is  shown  In  Figure  10.  This  grid  consists  of  114  points  In  the 
streamwlse  direction  and  50  points  In  the  normal  direction.  Grid  clustering 
In  the  streamwlse  direction  has  been  used  near  the  ogive-cylinder  and 
cyl 1 nder-boattal 1  junctions  as  well  as  the  base.  Figure  11  shows  the  expanded 
view  of  the  grid  In  the  base  region  and  again  shows  the  large  number  of  grid 
points  used  In  this  region. 

The  drag  components  (pressure  drag,  viscous  drag  and  base  drag)  for  this 
shape  are  shown  In  Figures  12  through  18,  Figure  12  shows  the  pressure  drag 
comparison.  At  higher  supersonic  Mach  numbers  .all  predictions  are  In 
reasonable  agreement.  The  pressure  drag  at  transonic  speeds  Is  under  pre¬ 
dicted  by  both  the  design  codes.  As  shown  In  Figure  13  the  pressure  drag 


13.  L.D.  Kays'er,  "Base  Pressure  Measurements  on  a  Projectile  Shape  at  Mach 
Numbers  for  0.91  to  1.20,"  US  Army  Ballistic  Research  Laboratory , 
Aberdeen  Proving  Ground,  Maryland,  ARBRL-MR-03353 ,  April  1884.  (AD 
A141341) 

14.  P.R.  Payne  and  R.M  Hartley,  "Afterbody  Drag,  Volume  1  -  Drag  of  Conical 
and  Circular  Arc  Afterbodies  without  Jet  Flow,"  Final  Report, 

DTNSRDC / ASED-80 /10 ,  Bethesda,  Maryland,  May  1980, 


obtained  from  the  Navier-Stokes  results  is  compared  with  experiment15  and  is 
In  excellent  agreement  with  experiment.  The  viscous  drag  is  compared  in 

Figure  14.  The  Navier-Stokes  result  is  In  good  agreement  with  the  MCDRAG 
prediction  at  all  Mach  numbers  from  .9  to  3.0.  The  NSWCAP  code  prediction  is 
In  reasonable  agreement  at  supersonic  Mach  numbers  and  Is  underpredicted  at 
all  transonic  Mach  numbers.  Comparison  of  the  base  drag  Is  shown  In  Figure 
15.  As  shown  In  this  figure,  the  Navier-Stokes  code  prediction  Is  In  good 
agreement  with  the  predictions  from  the  design  codes  at  higher  Mach  numbers  (M 
>  1.5).  Large  discrepancies  are  seen  between  the  predictions  at  transonic 
speeds.  Both  design  codes  MCDRAG  and  NSWCAP  grossly  overpredict  the  base  drag 
In  this  speed  regime.  The  Navier-Stokes  result  Is  compared  with  experiment13 
and  the  MCDRAG  result  In  Figure  16.  Although  there  Is  a  small  discrepancy 

between  the  numerical  result  and  the  experimental  data,  the  trend  l.e. , 

decrease  and  Increase  of  base  drag  with  Mach  number  seen  experimentally  Is 

clearly  observed  In  the  computational  results.  MCDRAG,  on  the  other  hand, 
does  not  show  the  correct  trend  and  grossly  overpredicts  the  base  drag  except 
at  H  *  1.2.  NSWCAP,  predicts  even  higher  base  drag  at  transonic  speeds  and 
thus,  Is  In  worse  agreement  with  the  Navier-Stokes  result  and  the  experiment. 
The  computational  result  and  the  experiment  show  a  negative  base  drag  (or 
thrust)  at  Mach  numbers  0.94,  0.96  and  0.98  which  Is  not  predicted  by  the 
design  codes. 

Figure  17  shows  the  variation  of  base  drag  with  Mach  number  for  both  SOC 
and  SOCBT  configurations.  These  are  the  computational  results  and  the  reduc¬ 
tion  In  base  drag  due  to  the  boattall  Is  clearly  seen.  Figure  18  shows  the 
total  drag  as  a  function  of  Mach  number.  The  Navier-Stokes  result  Is  compared 
with  the  design  code  predictions.  The  agreement  Is  good  at  all  Mach  numbers 
except  between  1.1  and  1.2  where  about  15-20%  discrepancy  Is  found  between  the 
design  code  predictions  and  the  Navier-Stokes  result. 

C.  M549  Projectile 

The  M549  Is  a  modern,  low  drag.  Army  artillery  shell.  The  geometry  of 
this  shell  Is  shown  In  Figure  19.  It  has  a  3  caliber  ogive  nose,  an  approxi¬ 
mately  2  caliber  cylindrical  section,  and  a  .59  caliber,  7  1/2°  boattall. 
Certain  simplifications  have  been  made  on  this  shape.  The  flat  nose  was 
modeled  as  a  hemisphere  cap  and  the  rotating  band  was  eliminated.  As  a 
result,  a  modified  confl guratlon,  shown  In  Figure  20,  was  obtained  and  used 
for  the  numerical  computations. 

An  expanded  view  of  the  computational  grid  near  the  projectile  Is  shown 
In  Figure  21.  This  grid  has  114  points  In  the  streamwlse  direction  and  40 
points  In  the  normal  direction.  Grid  points  are  clustered  near  the  surface 
and  also  In  the  base  region  where  large  changes  In  the  flow  variables  are 
expected.  Figure  22  shows  the  velocity  vectors  In  the  base  region  obtained 
from  the  Navier-Stokes  calculations.  The  reel rculatory  base  flow  Is  clearly 
evl dent. 


15.  L.D.  Kayeer,  "Surface  Pressure  Measurements  cm  a  Boattailed  Projectile 
Shape  at  Transonic  Speeds,"  US  Army  Ballistic  Research  Laboratory , 
Aberdeen  Proving  Ground ,  Maryland ,  ARBRL-MR-03 16 1 ,  March  1982.  (AD 
A113S20) 
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The  Individual  drag  components  (pressure  drag,  viscous  drag,  base  drag) 
and  the  total  drag  are  shown  In  Figures  23,  24,  25  and  26,  respectively.  The 
drag  components  obtained  from  Navler-Stokes  solutions  are  compared  with  the 
NSWCAP  and  MCDRAG  predictions.  The  pressure  drag  predicted  by  the  Navier- 
Stokes  code,  NSWCAP  code  and  MCDRAG  Is  In  good  agreement  at  low  Mach  numbers 
(.9  <  M  <  1.0).  MCDRAG  predicts  slightly  higher  pressure  drag  at  higher  Mach 
numbers  compared  to  the  Navler-Stokes  and  NSWCAP  predictions.  The  viscous 
drag  predicted  by  all  codes  Is  In  good  agreement.  As  for  the  base  drag,  the 
Navler-Stokes  result  agrees  very  well  with  the  NSWCAP  prediction  whereas 
MCDRAG  underpredicts  the  base  drag  by  as  much  as  20  to  25%.  The  total  drag 
comparison  Is  shown  In  Figure  26.  Again,  the  Navler-Stokes  result  Is  compared 
with  predictions  from  NSWCAP  and  MCDRAG  codes.  In  addition,  these  predictions 
are  compared  with  the  LCWSL16  data  base  and  are  all  In  good  agreement. 


IV.  CONCLUDING  REMARKS 

A  thin-layer  Navler-Stokes  code  has  been  used  to  compute  the  flow  field 
over  projectiles  Including  the  base  region.  Numerical  computations  have  been 
made  for  three  projectiles  (SOC,  SOCBT  and  M549)  for  various  Mach  numbers  and 
a  •<  0.  The  Individual  drag  components  (pressure  drag,  viscous  drag,  and  base 
drag)  and  total  drag  were  obtained  from  computed  solutions.  In  addition,  two 
design  codes  (NSWCAP  and  MCDRAG)  that  employ  semi -empl rl cal  techniques  were 
used  to  predict  the  drag  for  projectiles.  Predictions  from  the  Navler-Stokes 
code  and  the  design  codes  were  compared  with  experiment  and/or  data  base 
results  where  available. 

The  viscous  drag  predicted  by  the  design  codes  Is  generally  In  good 
agreement  with  the  prediction  from  the  Navler-Stokes  code.  However,  signifi¬ 
cant  discrepancies  In  the  prediction  of  pressure  drag  and  base  d*-ag  are 
found.  When  pressure  drag  Is  underpredicted,  the  base  drag  Is  overpredicted 
and  vice  versa.  Thus,  when  the  Individual  drag  components  are  added  up  to 
obtain  the  total  drag,  the  discrepancy  Is  usually  small  and  reasonably  good 
agreement  Is  found  between  the  predictions  of  total  drag  from  the  Navler- 
Stokes  code  and  the  design  codes. 


16,  A.  Loeb,  "Private  Cormunioation,"  Large  Caliber  Weapons  Systems 
Laboratory ,  ABDC ,  AMCCOM,  Dover,  Neu  Jersey. 
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Figure  4.  Variation  of  Pressure  Drag  Coefficient  with  Mach  Number 


Figure  5.  Variation  of  Viscous  Hrag  Coefficient  witn  Kach  Number, 
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Figure  7.  Rase  Drag  vs  Mach  Number  at  Transonic  Speeds, 


Figure  8.  Variation  of  Total  Oraq  Coefficient  with  Mach  Number, 
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at  ion  of  Pressure  Drag  Coefficient  with  Mach  Humber 
0,  SOCBT,  M  =  0.8  to  3.0 
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Figure  16.  Variation  of  Base  Drag  Coefficient  with  Mach  Number 
a  =  0,  SOCBT,  M  =  0.9  to  1.2 


Figure  17.  Base 
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Figure  19.  M549  Projectile 


Figure  20.  Computational  Model 
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Figure  26.  Total  Drag  Coefficient  vs  Mach  Number,  a  =  0,  M549 
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This  Laboratory  undertakes  u  continuing  effort  to  improve  t-V.'^flual ity  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
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1 .  BRL  Report  Number _ 

2.  Date  Report  Received 
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